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The traditional approach to accelerator optics, based mainly on clas- 
sical mechanics, is working excellently from the practical point of 
view. However, from the point of view of curiosity, as well as with 
a view to explore quantitatively the consequences of possible small 
quantum corrections to the classical theory, a quantum mechanical 
formalism of accelerator optics for the Dirac particle is being devel- 
oped recently. Here, the essential features of such a quantum beam 
optical formalism for any spin-i particle are reviewed. It seems that 
the quantum corrections, particularly those due to the Heiscnberg 
uncertainty, could be important in understanding the nonlinear dy- 
namics of accelerator beams. 
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1 Introduction 



Why should one bother about a quantum mechanical treatment of accelerator 
optics when the classical treatment works so well? This is a natural question 
indeed. There is no prima facie reason at all to believe that a quantum me- 
chanical treatment would be necessary to deal with any aspect of accelerator 
optics design. As has been rightly pointed out 0, primary effects in conven- 
tional accelerators are essentially classical since the de Broglie wavelength of the 
high energy beam particle is much too small compared to the typical apertures 
and the energy radiated is typically low and of long wavelength. However, it is 
being slowly recognized that quantum effects are still important due to demands 
on high precision in accelerator performance and ever increasing demands for 
higher beam energy, luminosity and brightness (jjj. Also, there is a growing 
feeling now that a complete picture of spin polarization can only be obtained 
on the basis of coupled spin <-> orbit phase-space transport equations (with and 
without radiation) and to include all the subtleties of radiation one has to begin 
with quantum mechanics since classical white noise models simply would not 
suffice for all situations Q . 1 like to add: 

• After all, accelerator beam is a quantum mechanical system and one may 
be curious to know how its classical behavior can be understood, in detail, 
from a quantum mechanical formalism based on the appropriate relativis- 
tic wave equation. 

• As has been revealed by recent studies jj) the passage from quantum the- 
ory to classical theory is not a simple affair, particularly when the system 
is a complicated nonlinear dynamical system with regular and chaotic re- 
gions in its phase-space. Since accelerator beams are such systems j| 
it is time that the quantum mechanics of accelerator optics is looked at 
seriously. 

Essentially from the point of view of curiosity, the axially symmetric mag- 
netic lens was first studied || based completely on the Dirac equation. Later 
works led to further insights into the quantum mechanics of spin-i particle 
beam optics. Quantum mechanics of the Klein-Gordon (spin-0) and nonrela- 
tivistic Schrodinger charged-particle beams were also studied m || . These works 
dealt essentially with aspects of ion optics and electron optical imaging (for an 
excellent survey of scalar electron wave optics see the third volume of the en- 
cyclopaedic three- volume text book of Hawkes & Kasper || ; this contains also 
references to earlier works on the use of the Dirac equation in electron wave 
optics problems like diffraction, to take into account the spinor nature of the 
electron). 

In the context of accelerator physics also, like in the case of electron and 
ion beam optical device technologies, the practice of design of beam optical 
elements is based mainly on classical physics. As is well known, various aspects 
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of accelerator beam dynamics like orbital motion, spin evolution and beam 
polarization, radiation and quantum fluctuations of trajectories, are analyzed 
piecewise using classical, semiclassical, or quantum theories, or a mixture of 
them, depending on the situation treated. Quantum mechanical implications 
for low energy polarized (anti)proton beams in a spin-splitter device, using the 
transverse Stern-Gerlach (SG) kicks, have been analyzed Jl(| on the basis of 
nonrclativistic Schrodinger equation. 

To obtain the coupled spin <-> phase-space transport equations, for the spin- A 
particle, one needs an appropriate quantum Hamiltonian. Such a Hamiltonian 
was stated, as following from the systematic Foldy-Wouthuysen (FW) trans- 



formation technique JjJ, by Derbenev and Kondratenko 12 in 1973 as the 



starting point of their radiation calculations but no explicit construction was 
given (such a Hamiltonian can also be justified |l3| using the Pauli reduction 
of the Dirac theory). The Derbenev-Kondratenko (DK) Hamiltonian has been 
used pi to construct a completely classical approach to beam optics, including 
spin components as classical variables. Now, a detailed derivation of the DK 
Hamiltonian has been given ^ and a completely quantum mechanical formalism 
is being developed 0] in terms of the 'machine coordinates' and 'observables', 
suitable for treating the physics of spin-^ polarized beams from the point of 
view of machine design. 

Independent of the DK formalism, recently p5[ fl6| we have made a beginning 
in the application of the formalism of the Dirac spinor beam optics, developed 
earlier (H-0) mostly with reference to electron microscopy, to accelerator optics 
to understand in a unified way the orbital motion, SG kicks, and the Thomas- 
Bargmann-Michel-Telegdi (TBMT) spin evolution. Here, I present the essential 
features of our work, done so far, on the quantum beam optical approach to 
accelerator optics of spin-i particles based on the Dirac equation. 



2 Quantum beam optics of the Dirac particle 

Our formalism of quantum beam optics of the Dirac particle, in the context of ac- 
celerators, is only in the beginning stages of development. So, naturally there are 
several simplifying assumptions: We deal only with the single particle dynam- 
ics, based on the single particle interpretation of the Dirac equation, ignoring all 
the inter-particle interactions and statistical aspects. Only monoencrgetic beam 
is considered. The treatment is at the level of paraxial approximation, so far, 
though the general framework of the theory is such that extension to the case 
of nonparaxial (nonlinear) systems is straightforward. Only time-independent 
magnetic optical elements with straight axis are considered. Electromagnetic 
field is treated as classical. And, radiation is ignored. 

Thus, we are dealing with elastic scattering of the particles of a monoen- 
ergetic beam by an optical element with a straight axis along the z-direction 
and comprising a static magnetic field B = curl A. Hence, the 4-component 
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spinor wavefunction of the beam particle can be assumed to be of the form 
\l/(r, t) = i/j(r) ex-p(—iEt/h), where E is the total (positive) energy of the parti- 
cle of mass to, charge q, and anomalous magnetic moment /i a . The spatial part 
of the wavefunction ip(r) has to obey the time-independent Dirac equation 

Htp(r) = Eijj(r) , (1) 

where the Hamiltonian H, including the Pauli term, is given by 

H = (3mc 2 + cot ■ 7T - ^ a /3S • B , 
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with TV = p — qA — — i?iV — qA. Let p be the design momentum of the beam 
particle along the +z-direction so that E = + y / m 2 c 4: + c 2 p 2 . The beam is 
assumed to be paraxial: |pj_| <C \p\ — p and p z > 0. 

Since we are interested in studying the propagation of the beam along the 
+z-direction we would like to rewrite Eq. (Q) as 

d 

ih—tl;(r±;z)=Htp(r±;z). (3) 

So, we multiply Eq. (|l|) from left by a z /c and rearrange the terms to get 
H = -pf3xa z - qA z i + a z a ± ■ ir ± + (n a /c)/3a z '£ B, 



( & \ , / 1 0\ , E + mc 2 

Note that the matrix Px a z, coefficient of — p in TC, is diagonalized as follows: 
MiPxa^Ar 1 = p , M=^=(t + X a z ). (5) 



Hence, let us dchnc 



rjj = Mtp . (6) 



This turns Eq. (g) into 

— 

dz 



ih^-ip' = H'f// , W = MUM- 1 = - P p + £ + 0, (7) 
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with the 'even' operator £ and the 'odd' operator O given, respectively, by the 
diagonal and off-diagonal parts of 



£ + 



[-qA t l-Qi a /2c) 
x{(Z + Z- 1 )v±-B J 
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The significance of the transformation in Eq. (|J) is that for a paraxial beam 
propagating in the forward (+2) direction tp' is such that its lower pair of com- 
ponents are very small compared to the upper pair of components, exactly like 
for a positive energy nonrelativistic Dirac spinor. Note the perfect analogy: 
nonrelativistic case — > paraxial, positive energy — > forward propagation, and 
mc 2 — ► — p in Eq. (||) and Eq. (Q) respectively. 

Let us recall that the FW transformation technique jy] [TtJ is the most 
systematic way of analyzing the standard Dirac equation as a sum of the non- 
relativistic part and a series of relativistic correction terms. So, the application 
of an analogous technique to Eq. (Q) should help us analyze it as a sum of the 
paraxial part and a series of nonparaxial correction terms. To this end, we 
define the FW-likc transformation 



ipl = Slip' , Si = exp (-P0/2p) 



(9) 



The resulting equation for tpi is 



*&-5-^i= Wi» 'Hx = Si'H'Si 1 -ihSx—{Si l \ =-p/3 + 8 1 + 1 , 
dz oz 1 J 



£1 = £ /30 2 

2p 



(10) 



A series of such transformations successively with the same type of recipe as 
in Eq. (j^) eliminates the odd parts from TL' up to any desired order in 1/p. 
It should also be mentioned that these FW-like transformations preserve the 
property of tp' that its upper pair of components are large compared to the 
lower pair of components. We shall stop with the above first step which would 
correspond to the paraxial, or the first order, approximation. 

Since the lower pair of components of ip\ are almost vanishing compared to 
its upper pair of components and the odd part of Hi is negligible compared to 
its even part, up to the first order approximation we are considering, we can 
effectively introduce a two-component spinor formalism based on the represen- 
tation of Eq. (|To|). Naming the two-component spinor comprising the upper 
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pair of components of ipi as ip" and calling the 2x2 11-block element of Hi as 
H" it is clear from Eq. ( |l0| ) that we can write 

oz 

H" « (- P -qA z + ^-Ttl)--{(q + e)B z S z + 7 eB x -S x } , 
V 2p J p 

with tt\ = fr x + 7r^ , e = 2m/j„/?i, r y~E/mc 2 , S^hcr/2. (11) 

Up to now, all observables, field components, etc., are defined with reference 
to the laboratory frame. But, as is usual in accelerator physics, we have to 
define spin with reference to the instantaneous rest frame of the particle while 
keeping the other observables, field components, etc., defined with reference to 
the laboratory frame. For this, we have to transform ip" further to an 'accel- 
erator optics representation', say, ipA = Tai/j" ■ The choice of Ta is dictated by 
the following consideration. Let the operator O correspond to an observable O 
in the Dirac representation (Eq. (pj) or Eq. @). The operator corresponding to 
O in the representation of Eq. < \l l| ) can be taken to be given by 

6" = the hcrmitian part of the 11-block element of (SiMOM^S^J . 

(12) 

The corresponding operator in the accelerator optics representation will be 

6 A = the hermitian part of (TaO"T a ~ 1 \ . (13) 



The operator 

/ _ e 2 (cr-7r7r+7rcr.7r) c7T \ 

S {R) =~h\ 2E(E+m<?) E_ / 14 n 

° o c7r _ . e 2 (cr.7r7r+7rcr.7r) I ' \ L *> 

\ — _<TH 2E(E+mc>) / 

corresponds to the rest-frame spin in the Dirac representation [ fl8[ . We demand 
that the components of the rest-frame spin operator in the accelerator optics 
representation be simply the Pauli spin matrices, i.e., Sa ~ 7wr/2, up to the 
first order (paraxial) approximation. This demand leads to the choice 

ipA = T A ip" , Ta = cxp {-i {fc x cry - n v a x ) /2p} . (15) 

Now, finally, with the transformation given by Eq. ( |l5| ) , the desired basic equa- 
tion of the quantum beam optics of the Dirac particle becomes, up to paraxial 
approximation, 

ih-^-ipA = Ha^a , Ha ~ ( -p - qA z + ^-Tr±] + — £t ■ S , 
oz \ 2p ) p 

with CI = — \qB + e (Bu + jB±) } . (16) 
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Here, Bu and B± are the components of B in the +z-direction (the predominant 
direction of motion of the beam particles) and perpendicular to it, unlike in 
the usual TBMT vector f2 in which the components B\\ and B± arc defined 
with respect to the direction of the instantaneous velocity of the particle. The 
quantum beam optical Hamiltonian Ha is the beam optical version of the DK 
Hamiltonian in the paraxial approximation. To get the higher order corrections, 
in terms of tt±/p and h, we have to go beyond the first FW-like transformation. 

It must be noted that while the exact quantum beam optical Dirac Hamilto- 
nian H (see Eq. (||)) is nonhermitian the nonunitary FW-like transformation has 
projected out a hermitian Ha (see Eq. (|l6|)). Thus, for the particle that survives 
the transport through the system, without getting scattered far away, ipA has 
unitary evolution along the z-axis. Hence, we can normalize the two-component 
ipA> a t an y z , as (tpA(z) \iPa(z)) = J f dxdyip A ipA = 1- This normalization will 
be conserved along the optic (z) axis. Then, for any observable O represented 
by a hermitian operator Oa, in the accelerator optics representation (Eq. (|l6|)), 
we can define the average at the transverse plane at any z as 

(0)(z) = (^a(z)\OaHa(z)) =J J dxdy^ A 6 A ipA- (17) 

Now, studying the z-evolution of (0}(z) is straightforward. Integration of 
Eq. ( |l6| ) gives 

il >A {z') = U{z l ,z)^ A {z), (18) 

where the unitary z-evolution operator U(z' , z) can be obtained by the standard 
quantum mechanical methods. Thus, the relations 

(0)(z') = (^a(z)\ U(z', zyd A U(z', z) \^ A (z)) , (19) 

for the relevant set of observables, give the transfer maps for the quantum 
averages (or their classical values a la Ehrenfest) from the plane at z to the 
plane at z' . In the classical limit this relation (Eq. ( [T9|)) becomes the basis 
for the Lie algebraic approach to classical beam optics |19|. The Lie algebraic 
approach has been studied |]20f in the context of spin transfer map also using 
the classical formalism. 

The main problem of accelerator optics is to know the transfer maps for 
the quantum averages of the components of position, momentum, and the rest- 
frame spin between transverse planes containing the optical elements. We have 
already seen that the rest-frame spin is represented in the accelerator optics 
representation by the Pauli spin matrices. Let us take that the observed posi- 
tion of the Dirac particle corresponds to the mean position operator of the FW 
theory |ll| or what is same as the Newton- Wigner position operator |21j| . Then, 
one can show, using Eq. ( |T2| ) and Eq. (|l|), that in the accelerator optics repre- 
sentation the transverse position operator is given by r±_ up to the first order 
approximation (details are given elsewhere p2[). For the transverse momentum 
in free space the operator is p ± in the accelerator optics representation. 
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3 An example: the normal magnetic quadru- 
pole lens 

For an ideal normal magnetic quadrupole lens of length L comprising the field 
B = (Gy,Gx,0), corresponding to A = (0,0, \G (y 2 — a; 2 )), and bounded by 
transverse planes at Zj and Zf = Zi + L, the quantum accelerator optical Hamil- 
tonian (Eq. ([l6|)) becomes 

-p + -^p\ , for z < Zi and z > Zf , 
n A = { -p+ fqG (x 2 - J/ 2 ) + l p pl {Wx + xay) t (20) 

for Zi < z < Zf . 

Note that the choice of A in a different gauge will not affect the average values 
defined by Eq. @. Now, using the formalism of the previous section, it is 
straightforward to find the desired transfer maps in this case (details are found 
elsewhere JT|, H))- The results are: with V = (l + ie)GLh/2p 2 , K = ^JqG/p, 
A = h/p, ()i = ( )(zi), and ( )/ = ( )(z f ), 

(x)f \ ^ ( cos KL sin KL 



(Px) f J \ —pK sin KL cos KL 

(x)i \,f (cos KL - l){a y ) i /K 2 L 
fe>, J V \ -(psmKL){a y )jKL 
(V)f \ J cosh KL ^sinhKL 
(Py)f J ~ \ pKs'mhKL cosh KL 

(V)i \ , ( ~{coshKL-l){a x )jK 2 L 



(Pv)i J V -(psmhKL)(a x ) l /KL 
Atttj J / sin if L \ . / sinh if L \ 



if 2 Lp J**"* y x 2 Lp 

Obviously we have obtained the well known classical transfer maps (matri- 
ces) for the transverse phase-space coordinates, and more, i.e., the transverse 
SG kicks jl(| ^3|. The longitudinal SG kick, which has been proposed as a 
better alternative to the transverse SG kicks for making a spin-splitter device 
to produce polarized (anti)proton beams, can also be understood Jl5| , p2| using 
the present quantum beam optical formalism. The skew magnetic quadrupole 
can also be analyzed p3 in the same way as here. 
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4 Conclusion 



In summary, we have seen how one can obtain a fully quantum mechanical for- 
malism of the accelerator beam optics for a spin-^ particle, with anomalous 
magnetic moment, starting ab initio from the Dirac-Pauli equation. This for- 
malism leads naturally to a unified picture of orbital and spin dynamics taking 
into account the effects of the Lorentz force, the SG force and the TBMT equa- 
tion for spin evolution. Only the lowest order (paraxial) approximation has 
been considered in some detail, with an example. It is clear from the general 
theory, presented briefly here, that the approach is suitable for handling any 
magnetic optical element with straight axis and computations can be carried 
out to any order of accuracy by easily extending the order of approximation. 
It should be emphasized that the present formalism is valid for all values of 
design momentum p from the nonrelativistic case to the ultrarelativistic case. 
The approximation scheme is based only on the fact that for a beam, consti- 
tuted by particles moving predominantly in one direction, the transverse kinetic 
momentum is very small compared to the longitudinal kinetic momentum. 

We hope to address elsewhere |2^] some of the issues related to the con- 
struction of a more general theory overcoming the limitations of the present 
formalism. With reference to the inclusion of multiparticle dynamics within the 
present formalism, it might be profitable to be guided by the so-called thermal 
wave model which has been extensively developed [ p5| in recent years to account 
for the classical collective behavior of a charged-particle beam, by associating 
with the classical beam a quantum- like wavefunction obeying a Schrodinger-Zifce 
equation with the role of h played by the beam emittance e. 

To conclude, let me emphasize the significance of the quantum formalism 
for beam optics. The following question has been raised: What are the uses 
of quantum formalisms in beam physics? jlj Of course, as we have seen above, 
we understand the quantum mechanics underlying the observed classical be- 
havior of the beam: when the ^-dependent quantum corrections are worked out 
through the higher order FW-like transformations it is found that they are really 
negligible. In my opinion, quantum formalism of beam optics has more signif- 
icant uses. To see this, let me cite the following cases: (1) Recently there is a 
renewed interest 1 26 in the form of the force experienced by a spinning relativis- 
tic particle in an external electromagnetic field. Such studies are particularly 
important in evaluating the possible mechanisms of spin-splitter devices |p3f . 
A thorough analysis based on general Poincare covariance, up to first order in 
spin, shows that classical spin-orbit systems can be characterized p7j , at best, 
by five phenomcnological parameters. This just points to the fact that spin be- 
ing essentially a quantum aspect one must eventually have a quantum formalism 
to understand really the high energy polarized accelerator beams. (2) A look 
at Eq. ([It]) shows that the form of the transfer map for the quantum averages 
of observables will differ from the form of the corresponding classical maps by 
terms of the type (f(Oi, O2, ...,))— /((Oi), (O2), ■ ■ • , ) 7^ 0; the differences are 
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essentially due to the quantum uncertainties associated with ipA a t the initial z%. 
For example, a term like x 3 in a classical map will become in the corresponding 
quantum map (x 3 ) = (x) 3 + 3(x) ((x — (x)) 2 ) + ((x — (x)) 3 ) which need not van- 
ish even on the axis (where (x) — 0). It is thus clear that, essentially, quantum 
mechanics modifies the coefficients of the various linear and nonlinear terms in 
the classical map making them dependent on the quantum uncertainties asso- 
ciated with the wavefunction of the input beam; actually, even terms absent in 
the classical map will be generated in this manner with coefficients dependent 
on the quantum uncertainties as a result of modifications of the other terms (a 
related idea of fuzzy classical mechanics, or fuzzy quantum evolution equations, 
occurs in a different context [p8| ). This quantum effect could be significant in 
the nonlinear dynamics of accelerator beams. This effect is relevant for spin dy- 
namics too. For example, it is seen in Eq. ( pl| ) that the spin transfer map is not 
linear in its components, in principle, even in the lowest order approximation, 
since terms of the type, say, (xS z ) i , (p x S z ) v etc., are not the same as (x)i (Sz)^ 
(Px)i(Sz) v etc., respectively, in general. 
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